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ABSTRACT 

A generalized approximate equation for duct lining sound attenuation is 
presented. The specification of two parameters, the maximum possible attenua- 
tion and the optimum wall acoustic impedance is shown to completely determine 
the sound attenuation for any acoustic mode at any selected wall impedance. The 
equation is based on the nearly circular shape of the constant attenuation contours 
in the wall acoustic impedance plane. For impedances far from the optimum, the 
equation reduces to Morse* s approximate expression. The equation can be used 
for initial acoustic liner design. Not least important is the illustrative nature of 
the solutions which provide an understanding of the duct propagation problem 
usually obscured in the exact calculations. Sample calculations using the approxi- 
mate attenuation equation show that the peak and the bandwidth of the sound atten- 
uation spectrum can be represented by quite simple functions of the ratio of actual 
wall acoustic resistance to optimum resistance. 



ATTENUATION OF SOUND IN DUCTS WITH ACOUSTIC TREATMENT - 


A GENERALIZED APPROXIMATE EQUATION 
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SUMM.\RY 

A generalized approximate equation for duct lining sound attenuation is pre- 
sented. The specification of two parameters, the maximum possible attenuation 
and the optimum wall acoustic impedance is shown to completely determine the 
sound attenuation for any acoustic mode at any selected wall impedance. The 
equation is based on the nearly circular shape of the constant attenuation con- 
tours in the wall acoustic impedance plane. For impedances far from the opti- 
mum, the equation reduces to Morse’ s approximate ejq)ression. 

The equation can be used for initial acoustic liner design. Not least impor- 
tant is the illustrative nature of the solutions which provide an imderstanding of 
the duct propagation problem usually obscured in the exact calculations. 

Sample calculations using the approximate attenuation equation show that 
the peak and the bandwidth of the soimd attenuation spectrum can be represented 
by quite simple functions of the ratio of actual wall acoustic resistance to opti- 
mum resistance. 


INTRODUCTION 

The approximate sound attenuation equation devele>ed by Morse (ref. 1) 
which is valid for nearly hard walls has been extremely useful and produces 
quite good results over its range of validity. This equation is easy to use re- 
quiring only liner resistance, reactance, and duct length to diameter ratio as 
inputs. Unfortunately, as liner resistance is reduced and the exact calculations 
using the wave equation predict higher sound attenuaticms, the approximate equa- 
tion becomes inaccurate. 

The purpose of this paper is to generate an approximate sound attenuation 
equation which is as simple as possible in form but which will adequately com- 
pare with some exact calculations in the lower resistance or high attenuation 
region of tiie wall impedance plane. The approximate equation is anchored on 
the optimum impedance and the maximum possible attenuation associated with 


I 


2 

this impedance. This optimum has been extensively studied and reported for 
single (ref. 2) and multiple modes (ref. 3) without flow, multiple modes with 
flow (ref. 4), for spinning modes with uniform flow (ref. 5), and with boundary 
layers (ref. 6) . The mathematical implications of, and a technique for deter- 
mining this optimum impedance have been reported (refs. 7 and 8). A complete 
correlation of the optimum point for arbitrary spinning modes in a circular duct 
with uniform flow is included in this paper. 

With this well established optimum point as an input to the approximate at- 
tenuation equation, the liner attenuation at any other wall impedance can be cal- 
culated. The model is used to generate illustrative examples which provide in- 
sight into the operation of acoustic suppressors. Design tools or "rules of 
thumbs* are also generated for the peak spectral attenuation and the attenuation 
bandwidth. 


SYMBOLS 

Bg optimum resistance coefficient 

B^ optimum reactance coefficient 

b resonator backing depth, m 

C speed of sound, m/sec 

D circular duct diameter, m 

AdB sound power attenuation, dB 

AdB^ maximum possible sound power attenuation which is a function of 
frequency and other variables, dB 

AdB ^ P value of at the frequency of peak sound attenuation, dB 

AdBp peak value of sound power attenuation 

F function of maximum possible attenuation (see eq. (8)) 

f frequency, Hz 

frequency of peak spectral attenuation, Hz 
f^ upper frequency at which half of peak attenviatian is attained, Hz 

vq>per frequency at which half of peak attenuation is attained, Hz 

i 

acoustically treated duct length, m 
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Mq uniform steady flow Mach number in duct 

m spinning mode lobe pattern number 

Q 1 + iMQ((T + ir) 

R amplitude of eigenvalue o; 

(ft radius of constant attenuation contour in wall impedance plane 

Tq circular duct radius, m 

a complex eigenvalue (a = Re^*^) 

a_ complex eigenvalue for mode with m circumferential lobes of radial 

m,#i , 

order /n 

^ ratio of maximum possible to actual sound power attenuation at a par- 

ticular frequency 

/3p value of /3 at the frequency of peak attenuation 

y - 1 or /3 - 1 

5 boundary layer thickness, m 

T) frequency parameter, fD/C 

T) Tj at the frequency of peak attenuation 

6 specific acoustic resistance 

6 resistance coordinate of the center of an equal attenuation or damping 

c 

contour in the specific acoustic impedance plane 

optimiun specific acoustic resistance which is a function of frequency 
and other variables 

6 value of 6 at the frequency of peak sovmd attenviation 

mp iti 

radial mode number 
I normalized frequency, f/f^ 

a attenuation coefficient 

maximum possible value of cr for given set of conditions 
T propagation coefficient 

(p phase of eigenvalue at 

X specific acoustic reactance 
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X- reactance coordinate of the center of an equal attenuation or damping 
c 

contour in the specific acoustic impedance plane 

optimum specific acoustic reactance which is a function of frequency 
and other variables 

value of Xjjj at the frequency of peak sound attenuation 
CO circular frequency, rad/sec 


DEVELOPMENT OF THE APPROXIMATE EQUATION 

The approximate equation that will be presented for the equal damping con- 
tours was suggested by the shape of these contours as shown in figure 1. The 
solid line contours represent exact calculations for a plane wave input into the 
acoustic liner (ref. })• Althou^ the contours are cpiite irregular for high attenu- 
ation they become more nearly circular as the attenuation is reduced. Even for 
hlg^ attenuations, a circular contour would usually be sufficient for first approxi- 
mations. NOkice that, as the attenuation is reduced, the center of the contours 
move to higher resistance. All contour plots seem to show this behavior. Re- 
sults for spinning modes in ducts with sheared flow will be presented later in 
this paper. Also shown on figure 1 are the results calculated from Morse* s ap> 
proximate equation which is valid for nearly hard walls or for very low attenua- 
tions (less than 10 dB attenuation for the conditions of fig. 1). If this system of 
circular contours (Morsels approximate results) could be shifted and to the 
left, perhaps an adequate approximation to the exact calculations could be 
achieved even near the optimum impedance. This Impedance shift of the circur 
lar contours is exactly what the approximate expressions of this p{^r are in- 
tended to accomplish. 


Approximate Equation of Morse 

The approximate sound attenuation esqpression, representing the dashed con- 
tours of figure 1, which is valid for a plane wave entering a circular duct with 
nearly hard walls (low attenuation) is given by (ref. 1) 

w ^ yP) (1) 

where AdB is the sound power attenuation in decibels, 6 and x the liner 
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specific acoustic resistance and reactance, L is the length of acoustic treat- 
ment, and D is the diameter of the circular duct. The no- flow (Mq = 0) condi- 
tion is represented by equation (1). Attenuation is understood here and the usual 
minus sign is omitted. Equation (1) can also be written 


which represents a system of circular contours with the centers located at 


, 8.7 L/D 
AdB 


( 3 ) 


Xc 


= 0 


and of radius 


AdB 


( 4 ) 


(5) 


Notice that the approximate expressions already have some of the properties 
noted in the exact calculations: as the attenuation is reduced the center of the 
circle moves to higjher resistance (9 J and the radius (A) increases. 

V 


Improved Approximate Equati(m 

When making the translation of the approximate equal damping contours in 
the impedance plane, certain restrictions are set tqpon the equations. 

(1) As the maximum possible attenuation is approached the contour radius 
must approach zero. 

(2) At very low attenuations the results should approach equation (1) which 
is known to be valid for low attenuations. 

(3) The contour radius must always be less than the resistance ordinate of 
the circle center (negative 6>b not allowed). 

A simple system of equations which satisfies the above restrictions is 

»c - V ^ Fy (6) 

$ y 

(R-«S^+Fy 
1 +y 


( 7 ) 
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where is the optimum resistance and where 

p _ 8.7 L/D 

AdB^d + M„)2 


( 8 ) 


and 


y = 


AdB^ 
rn 

AdB 


1 


(9) 


with the subscript m on ^dB^ implying the maximum possible attenuation for 
the given conditions. Equations (6) to (9) can then be iised in 


to provide a reasonable estimate to the exact calculations. The above results 
were obtained by fitting 0 and (R to the results obtained from figures 4 and 5 
of reference 3 with the Mach number correction determined from reference 5. 

Note that a plane wave incident upon a soft walled section may produce several 
modes in the suppressor. The energy balance between, and the effect (vx>n over- 
all sound attenuation of, these modes depend iqson the frequency, duct size (L and 
D), and the wall impedance. These approximate expressions (eqs. (6) to (10)) do 
not precisely fit the exact calculations over the wide range of conditions considered 
because of the presence of the several modes. In a later section a more precise 
expression which is very adequate for single modes will be developed. The be- 
havior of multimodal distributions may then be approximated by repeated use of 
the single mode equations for several different modes. However, much informa- 
tion has been gained by the use of the cruder e3q>ressi<ms given above and the de- 
velopment will proceed with these expressions. 

Returning to equations (6) to (10), it is interesting to lo<dc at their limits. 

Note that as AdB -* then y -► 0, 0^ -* and (R 0 as required. Also 

as AdB 0, y w AdB^^^/AdB and 


8.7VD 
AdB(l + Mq)^ 


( 11 ) 


thus 


t.iUn 

AdB<l + M,'* 




( 12 ) 
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<R 8J 

AdB(l + Mq)^ 


( 13 ) 


and equation (10) becomes 


^ 8. 7 L/D 

AdB(l + Mq)^ 

X 

which is like equation (2) (when Mq = 0) except for the Xj^ term. Thus the re- 
sult for very small attenuation is like Morse' s approximate expression except for 
the slight shift in reactance. 

Notice that all that is now required to cc mpletely specify the damping 
contours throughout the wall impedance plane (using eqs. (6) to (10)) is the opti- 
mum resistance (®jj,) and reactance (Xjjj) maximum possible attenuation 

(AdBjjj). Thus off-optimmn liner performance can be estimated using only the 
optimum properties and the liner impedance behavior. These optimum properties 
will be specified in a later section. 

The inverse problem is often of interest: given the resistance and reactance 
of an acoustic liner, what is the damping? Equations (6), (7), and (10) can be 
combined to give 


^ (X - x^)" « 


8. 7 L/D 


AdB(l + Mq)' 


(14) 


® ’ 20(0^ + Fy) +(x 


X + 




(1 +Y) 


2L [e iLLM + 2Fyl = 0 
m * 


(1 +Y) 


/ = 


(15) 


If the damping ratio is defined as 


in 

AdB 


= y + 1 


(16) 


then equation (15) is 

^ * 266m ■ 26FO - D + (x - X^)^ + ®m ^m ^ = 0 (17) 

Note tha-‘ to solve for attenuation (AdB), /3 must first be determined from equa- 
tion (17) which is cubic in 0. lastead of using a closed form (but cong>licated) 
cubic equation technique, a procedure taking increments in 0 was used (starting 
with ^ " 1) until equatitm (17) was satisfied. Ai^roadmate solutions to equation 
(16) or (17) can also be derived for small y or large 0, respectively. For 
example, when 0 is large, 
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+ 28F + 2e„F + (x - 
2 OF 


( 18 ) 


or 


AdB » 


17.4 g(L/D) 


. 200^ + 20F + 20^F + (X - X^)^ d + Mq)^ 
which provides a considerable improvement over equation (1). 


(19) 


INPUTS REQUIRED FOR THE APPROXIMATE ATTENUATION EQUATION 

In this section the inputs required for the approximate attenuation equations 
will be presented. These are intended mainly for the sinj^e mode equations de- 
veloped later in the paper. For the simpler equations of the previous section 
(plane wave) the inputs will be developed as they are used in the next section. 

Inspection of equations (8), (16), and (17) shows that M^, L, D, AdB^, 

0^, and Xjj^ must be specified before the approximate equations can be used 
to calculate AdB. Also the frequency is necessary to determine the frequency 
parameter 

Tj » ^ (20) 

c 


since 0^, and Xj^ are functions of t). Mq, L, D, and ij will be 

quite easily determined since they will be dictated ^ the mass flow, size, and 
frequency of the noise source and thu allowable duct lei^dl^. The results which 
follow have been determined by empirical correlations of the optimum resis- 
tance, reactaiice, and eigenvalues obtained trom exact but tedious calculations 
involving sin^e mode soluticms with lobe number from 1 to 20 and tqp to the 
tenth radial mode. The eigenvalues will be used to calculate maximum possible 
attenuaticm (AdB^j^). Circular ducts without centerbodles are considered here. 

The optimum resistance and reactance are determined from (ref. 5) 

*m ■ «« 

and 


where 
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Q = 1 + 1M(<t + It) 



(23) 


(24) 


with a and r representing the damping and propagation coefficients of the mode 
(part of the complex wave number) and the eigenvalue is 

a = (25) 

The optimum resistance and reactance coefficients (B^ and B^) will be pre- 
sented belm^ . Double subscripts (indicating specific modes) have been left off 
for brevity ,on 0^, X|jj» B^, Q, a, t, and a) but it must be recognized that 
we are considering a particular mode in the correlations given here. For in- 
stance, a should really be considered as a_ u where m is the circumferen- 

in, 

tlal mode number and p is the radial mode number. Using this more complete 
notation, the eigenvalue correlations can be e^q^iessed as 

Rm.p ” V 1 + 0* 076 m^f^l (26) 


<P^ * 

'^m.p 


<P 


m, 1* 


-(M- » ■ 
,2(Vp+!n/7) 


(27) 


where 


, « m + 2. 247 m^''^ + 1. 521 

nil i 


(28) 


and 


<P 


m, 1 


(m + 2)®' ® 


(29) 


EquatiOB (28) is good for m 0. If m > 0, use Rq 278. The resistance 
and reactance coefficients can be ‘estimated from, 
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Bg « — (30) 


B 


X 


« 0.92 

+ 3. 17 


(31) 


Note that equations (26) to (31) are valid only at the optimum impedance which is 
the only input point needed for the approximate attenuation equations. Finally 
the maximum possible attenuation is calculated from 


AdBj^= 17.4 7rrjtT^(L/D) 


(32) 


where again attenuation being understood, the minus sign has been dropped and 
the subscript on sigma signifies maximum possible attenuation coefficient. 

If equati<ms (25) to (29) are used with p » 1 to define the eignevalue, and 
equation (24) is used to calculate the damping coefficient then equation (32) 
will provide an approximate reproduction of figure 5 in reference 5 for the damp- 
ing of the least attenuated spinning mode. 

A correction may be required for the optimum resistance and reactance in 
inlets if the boundary layer thickness is significant. This correction is corre- 
lated in reference 6 fox' well wot- cm modes. By well cut- cm it is implied that 



For modes that are approaching cut-off, a boundary layer refraction correlation 
has not yet been developed. However as cut-ui is approached, the wave fronts 
travel more transverse to the velocity gradients and the refraction effects should 
be reduced. 

The correlations given in dxis section are mainly intended to be used with 
the improved single mode iqpproximate attenuation expressions developed later 
in this paper. In the next section where plane wave inputs are considered, even 
more simplified approximationn will be used. 
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SAMPLE CALCULATION - PEAK ATTENUATION AND BANDWIDTH 

Illustrative Example 

Before the calculations for peak attenuation and bandwidth are made, it is 
instructive to illustrate the frequency dependent behavior of the equal attenuation 
contours and of a typical liner impedance. Such an Illustration is shown in fig- 
ure 2 where the attenuation contours are shown for three frequencies each an 
octave apart. The properties of the optimum points shown here imply well cut- on 
behavior for all three frequencies. These properties evolve from equations (20) 
to (22) which show that 


OC f 

(34) 

a -f 

(35) 


and equations (24) and (32) which show that (see ref. 5) 

AdBm a 1/f (36) 

The OT/timiun points are shown as plus signs; and these, as well as the equal 
damping contours, are labeled with their relative attenuations. Superimposed 
upon the damping contours are some sample liner resistances and reactances. 
These liner impedance loci assume that the resistance (6) is constant with fre- 
quency while the reactance of a Helmholtz resonator at low frequracies behaves 
as 

(37) 

which implies that back cavity stiifnesf is controlling the reactance. 

First consider the lower resistance locus (open symbols) which pass through 
the optimum point at the center frequency. As frequency increases the liner 
impedance moves to the right toward less negative reactance (eq. (37)) addle the 
optimum point and its related contours move to the uppe.r left toward more posi- 
tive resistance (eq. (34)) and more negative reactance (eq. (3f>)). As long as the 
reactance is to the left of its aasociated optimum (remember that each liner im- 
pedance point and each contour set have a frequency issociated with them) the 
damping is increasing with frequency as ever higher damping contoura are being 
cut When the reactance moves to the right ci its associated contour set the 
damping will f«il off and when X ” damping will occur. An example 

will illustrate these points. At f ■ 1/2 fp the far left wall lnq>edance point 
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(0=1, X = "2) is associated with the contours centered at the lower ri^t and a 
relative damping of 0. 52 dB is obtained. At f = fp, the liner impedance is coin- 
cident with the optimum point of the center contour set and a relative damping of 
1 dB is obtained. For the far ri^t liner impedance (1, -0. 5) at f = use the 
contour set centered at the upper left to read 0. 19 dB. In this example three 
points have been established in the sound attenuaMon spectrum. More points 
could be similarly determined at other frequencies. 

Now if this exercise is repeated at a hi^er resistance of 0=3 shown by 
the solid symbols in figure 2 it is seen that the liner impedance never passes 
ttorou^ an optimum impedance point but tiie relative damping still has a relative 
peak at the center frequency where x ~ Table I reviews the above results 
with an intermediate resistance also included. 

Note that a modest increase in resistance (0 = 1,6) above the optimum 
causes the peak damping to be reduced only moderately while both the hi^ and 
low frequencies have improved attenuation. A further increase in resistance 
(0 = 3) improves the hi^er frequency damping but decreases the damping at both 
of the other frequencies. Further increases in resistance would result in falloff 
of attenuation at all frequencies. The trade-off between peak attenuation and 
bandwidth would be determined by the specific application being studied. 

An exact calculation using the wave equation and more exact liner proper- 
ties would be performing the same ritual as illustrated in figure 2 except that 
the significance of the visualization would be lost in the process. A more 
sophisticated approach would use many modes at a time meaning that multiple 
contour sets would be considered simultaneously. 


Equations and Approximations 

In this section the peak attenuation and the attenuation bandwidth with a 
plaiie wave input are determined to provide one illustration of the use of the ap- 
proximate attenuation equations. Both quantities are functions mainly of the 

ratio of actual resistance to the optimum resistance at the peak attenuation fre- 

2 

quency. The starting point is equation (17), vdiich when divided by 0j^^p, yields 


u 

0 




-]■ 


^ - Xm> ^ 

0?-_ 0L_ 0^_ ®mp 


- 1 ) 


, 2FW - 1) 


1 

L. m 


0 


mi 


= 0 (38) 


mp mp 'mp 

Recall that 0 is the suppressor resistance, 0^^ is the optimum resistance which 


i * 

i 


! 


i 



13 


is a function of frequency, and is the optimum resistance at the frequency 
of peak attenuation at which x ~ Xjjj- I® the development which follows (to eq. 

(43)), some very simplified interrelationships among the quantities in equation (38) 
will be used. The objective is to reduce equation (38) to a form which gives attenu- 
ation (throu^ P) as a function only of a frequency ratio, a resistance ratio, and a 
peak frequency parameter. Define the quantity 



(39) 


where fp is the frequency at which x = Xjjj- From equation (20) and (21) 


m 




^mp ’^p ^p 


(40) 


2 

where it was assumed that Q is not a function of frequency, which implies that 
either Mq = 0 (then Q = 1), or that the frequencies involved are all above cut-off 
leading to Q = l/(l + Mq) (see ref. 5). Let the reactance be controlled by the 
back cavity 

X « -cot — » — (41) 

c 27rltf 


where the second approximation is valid for small arguments oi the cotangent, 
and where b i. the liner backing depth. Also 


(42) 

2"btp f { 

where equation (41) was used to define x^p as x at 
and (22) yield 


B 0 0 

V =-_iLiS«->SL 

B. 3 


(43) 


where Bq « 3B^ is a roug^ estimate from table I of reference 6 for axisym- 
metrlc wave propagation. Then from equations (40), (42), and (43) 






2 


( 44 ) 
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If, instead of using the approximation in equation (41), the cotangent had been 
retained, the result would be 



(45) 


An approximation to the maximum possible attenuation results of reference 3 
can be expressed as. 


m 

L/D 


» ' 


65 


^1 + 38 


1/4 


Finally using equations (8), (21), and (46) 


7.5n 


1/4 


mp 


(46) 


(47) 


where « 1 and Q w i/(i + Mq) were used. For sufficiently hi^ frequencies 

(high i]p^) 


F 





i 

3 


(48) 


biserting equations (40) and (44) into (38) yields, 



+ 4 


\fmp/ ^ 


0 


(49) 


Equation (49) along with (47) or (48) (depending upon or frequency) can be 
used in two ways. ' Irst with $ » I the peak attenuation was calculated throue^ 
/3p as a function mainly of The inverse of ^p (AdBp/AdB^^ deter- 

mined from equation (49) Is plotted In figure 3. For liner resistance below the 
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optimum resistauce, a single curve is obtained for all frequency parameters 
(tJjj) and the result is 

AdBp = ^d.lB„p; (e<«mp) (5») 

mp 

a conveniently simple expression. Equation (49) can be seen to contain the fac- 
tor (/3 - Q/®mp) I ~ 1. which accounts for this simple relationship. 

For liner resistance above the optimum the peak attenuation is seen to be a 
weak function of frequency parameter (r]p). An adequate approximation is, 

AdBp^^AdB^p; (9>9„p) (51) 


These two relationships, equations (50) and (51), are extremely useful ’’rules 
of thumb” in preliminary suppressor design considerations and in interpretation 
of acoustic liner data. 

Some results of exact calculations obtained from reference 3 are included 
on figure 3 for comparison with the approximate results obtained from equa- 
tion (49). 

The second use of equation (49) was to fix p at some desired level and de- 
termine the frequency ratio (| = f/fp) at which this attenuation occurred. This 
was done to determine the attenuation bandwidth or the frequencies at which 
AdB = l/2 AdBp, for instance. This can be accomplished by the following. Let 

AdB » 1/2 AdBp (52) 

then 


AdB^ AdB„ AdB^„ AdB^ 

s SI B tap m_ 


AdB 

Use of equation (46) gives. 


AdB AdB. AdB. 


20 


AdB. 


m 


mp 


P AdB 


mp 


AdB. 


m 


! 1 + 


AdB 




mp 


38n 


1/4 


\l + 


« - (for large tj^$) 


(53) 


(54) 


The lower and iy)per frequencies (fj/fp «nd f 2 /fp) for which one-half peak 
attenuation occurs were calculated with both representations of the reactance 
(eq. (41)). The results are shown in figures 4 and 5 with all of the calculations 
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made with the cotangent representation of reactance except for the curves labeled 
• 'large rjp limit, lumped reactance." For large the cotangent representa- 
tion approaches the lumped reactance result. The breakaway from the ''lumped 
reactance" curve for > 1 in figure 5 occurs where the back cavity half 

wave frequency is being approached, the reactance is approaching infinity and 
the attenuation is going to zero. In a real liner where a hi^ damping resonator 
is used this cannot occur and the use of the lumped reactance curve is probably 
more nearly correct for constructing the attenuation spectrum. Some fall off in 
attenuation will be anticipated beyond this half wave point. 

Some interesting conclusions can be seen from figures 4 and 5. For liner 
resistance below the optimum < 1) bandwidth remains essentially constant 

since fi/f„ and fn/f-, do not change much. Recall from figure 3 that peak at- 
tenuation falls off as resistance is reduced below the optimum. Thus no im- 
provement in liner performance is gained from underdamping the liner. How- 
ever, for high resistance > 1)» althou^ peak attenuation is lost, band- 

width is gained since fj/fp is reduced and f 2 /fp is increased. Recall from the 
discussion of figure 2 that increasing the resistance above the optimum value to 
improve attenuation bandwidth can be pushed too far since at very hi^ resistance 
the entire attenuation spectrum will be reduced as resistance is increased. 

Some useful " rules of thumb" can also be determined for bandwidth. These 
are obtained from observation of the results in figures 4 and 5. For resistances 
below the optimum. 




For resistance above optimum 


Vp) 


opt 


mp 


opt 


(55) 

(56) 


(57) 

(58) 


The subscript opt implies the condition when 0 = The following correla- 

tions are offered at the optimum resistance for use with equations (55) to (58). 
These were obtained by plotting the results at 8/^u^p = 1 from figures 4 and 5 
against t)p. 


i 
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(59) 


(60) 

Althou^ the sample calculations made here used information about circular 
duct optima, the final results (peak attenuation and bandwidth "rules of thumb") 
woixld be expected to apply to other geometries as well. For other geometries 
only the constants (such as ^dB^) in the problem would change. The 

key behavior of the approximate attenuation contours would not change. This 
behavior includes the vmsymmetrical shape of the contour pattern, the movement 
of the contours in the impedance plane as frequency changes, and the fall- off in 
AdBjjj as frequency increases. One precaution must be obseiwed, however, and 
that is in the use of equations (59) and (60). The frequency parameter at peak at- 
tenuation (Tip) must be considered relative to the cut-off frequency parameter. 

For example when rjp = i, cut-off is approached in the example shown. For 
smaller values of t]^, approaches a constant (losing the l/f dependence). 

For the first radial mode of a 20- lobe spinning pattern in a circular duct, AdB^^^ 
versus r\ will flatten out below tj = 10 (see ref. 5). 

The three quantities calculated here (AdBp, f^, and fg) can be used to obtain 
a first approximation to the attenuation spectrum. Log AdB versus log f is 
used with straight lines connecting the three points to obtain a reasonable approx- 
imation to the attenuation spectrum. Of course, the entire attenuation spectrum 
could be calculated from the approximate attenuation equations by using a multi- 
tude of frequencies rather than just fj^, fp, and fg. 

These simplified notions about peak attenuation and bandwidth have been used 
to analyze the attenuation data for inlet suppressors with splitter rings mounted 
on full scale fans and engines. These results are reported in reference 9 and 
show fairly good agreement between theory and data. Splitter ring suppressors 
seem to be relatively Insensitive to spinning mode c<mtent (ref. 5) and can be 
adequately described using the least attomated radial mode, 

IMPROVED SINGLE MODE APPROXIMATE ATTENUATION EQUATION 

Due to the form of the approximate attenuation equations discussed in this 
paper, the attenuation is assured of being an accurate representation of the more 
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exact propagation calculations at the extremes of damping. At the optimum im- 
pedance the attenuation will become the maximum possible and far from the opti- 
mum Impedance the attenuation will approach zero. Between these extremes a 
good fit to exact calculations was obtained by the empiricism inserted in the ap- 
proximate equations. As discussed previously, a plane wave input to a soft 
walled duct leads to a multimodal solution in the soft duct. The balance among 
the modes in determining the overall attenuation depends upon the wall imped- 
ance, the frequency and the liner geometry. A good fit of the approximate equa- 
tions to the exact calculations cannot be made over all conditions. A better ap- 
proach is to use approximate equations for each mode and to build up the total 
solution from the component parts. An improved single mode solution to pro- 
vide the basis for such an approach will now be developed. 

The resistance coordinate for the circular contour center was not altered, 

( 6 ) 

but the radius has been changed to, 


(ft = 


m' 


(V + l)(y + 2) 


+ Fy = ^ ™ 




+ - 1 ) 


(61) 


In addition, the reactance coordinate for the center of the constant attenuation 
contours can be represented as^ 


^ _ Xm ^ X|n 

® y + 1 

The definition of F has not been altered and is repeated as. 
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Technically the definition of F should have been altered since the solution for 
nearly hard walls results in. 
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Where m is the spinning pattern lobe number and is the mode eigenvalue 

for hard walls. Empirically this refinement (using m )) was found to 

be unnecessary in any reasonable vicinity of the optimum impedance (attenuations 
down to l/50 of the maximum possible) and this term was omitted from the de- 
velopment. For extremely low damping equation (63) mi^t be used with equa- 
tion (10) (with Xjjj = 0) if desired. 

Use of equations (6), (61), and (62) in equation (10) (with replacing Xjjj 
in eq. (10)) yields. 


20 |0jj^ + F(^ - 1) 




/ 3 (^ + 1 ) 


-n + i) 


m + 1 ) 


+ 2F(/3 - 1) = 0 

J 


(64) 


For circular ducts, the inputs (0^^, Xjjji calculated from the 

results given in a previous section (INPUTS REQUIRED FOR THE APPROXIMATE 
ATTENUATION EQUATION) and the attenuation (throu^ j3) can then be calculated 
at any resistance (0) and reactance (x) using equation (64). This equation is now 
seen to be of fifth power in p and a stepping or iteration approach must be used. 
Usually the ^ denominator of x^ can be omitted without much loss of accuracy. 

The equal attenuation contours (expressed for convenience by the -attenuation 
coefficient a rather than AdB, see eq. (32) for relationship) of the improved 
single mode approximate equation are compared to the exact calculations for a 
seven lobe spinning wave in figures 6 and 7. Figure 6 illustrates the first radial 
while figure 7 is for the fifth radial mode. The agreement between the approxi- 
mate and exact calciilations is seen to be quite good in both cases. The first 
radial mode would be considered well cut- on while the fifth radial is approaching 
cut-off (the left hand side of eq. (33) is »0. 46 which is not much less than 1). 

The rotation of equal attenuatian contours which occurs in the wall inq>edance 
plane due to in equations (21) and (22) has been accounted for and the entire 
contour pattern rotates with the optimum point. The exact equal damping c<m> 
tours of figures 6 and 7 were calculated as in reference 6 using a fairly thin 
boundary layer. Note the exact contours are tear-drop shaped for the higher at- 
tenuations but rapidly approach circular shape as the danqplng of the contours Is 
reduced. The shift of the equal contour centers toward hic^er resistance as 
dangling decreases is much more obvious in figures 6 and 7 dian in figure l. 


\ 
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CONCLUDING REMARKS 

Many equal damping contour plots have been investigated beyond the three 
that have been presented in this paper. These include a wide range of spinning 
lobe numbers and radial mode munbers, both with and without steady flow bovmd- 
ary layers. AU of these plots seem to have the bohavior described in this paper. 
The contours are nearly circular with the centers of these circles moving to 
higher resistance as the attenuation of the contour is reduced. 

With the contour properties and Morse* s approximate equation as guides, 
the approximate attenuation equations were developed. Two versions of the 
approximate equations ware pre sented. The first version has some contour 
broadening due to the multimodal nature of the exact solution to which it was fit. 
The second version, which is the favored version for future use, was generated 
to fit single mooes only. As a result this version is a more accurate represen- 
tation of more exact calculation procedures. 

Although the approximate attenuation equation was presented only as an 
alternate to more time consuming exact calculation procedures, such approxi- 
mate solutions may represent the most sensible approach for acoustic liner 
design. Due to the unknovm nature of the modal input to be used and the possible 
errors in physically obtaining the desired (analytical) value of wall acoustic im- 
pedance, it is possible that exact wave equation solutions are not justified. 

Using the approximate atti^nuation eqtiation the following " rules of thumb" 
were observed. 

1, For liner resistance above the optimum, the peak attenuation is Inversely 
proportional to resistance and t^e attenuation bandwidth increases with resistance. 

2. For liner resistance below the optimum, the peak attenuation is reduced 
proportionally with the resistance but the attenuation bandwidth remains constant. 
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